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We present a microscopic examination for the itinerant-localized duality model which has 
been proposed to understand anomalous properties of strongly correlated systems like the heavy 
fermions by Kuramoto and Miyake, and also useful to describe the anomalous properties of the 
high-T" c cupurates. The action of the duality model consists of two components of dynamical 
degrees of freedom. One of them is coherent itinerant degree of fermion and the other is 
incoherent localized one describing mainly the spin which has not been explicitly considered 
so far. We show that the thermodynamic potential of the strongly interacting Hubbard model 
can be rearranged in the form of duality model on the basis of renormalized perturbation 
expansion of the Luttinger-Ward functional if the one-particle spectral weight exhibits triple 
peak structure. We also examine the incoherent degrees of freedom described as a "localized 
spin" and show on the basis of the pertubation expansion that there exists commensurate 
superexchange-type interaction among the "localized spins" . 
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§1. Introduction 

The Fermi liquid theory offers us a very powerful tool for describing the low energy properties 
of strongly correlated fermions. While it was first devised by LandauEJ to analyse the physics 
of liquid 3 He, its general and profound value has been recognized through the application to the 
many areas of correlated fermion systems such as Kondo problem for instance. But it is also 
true that there are many phenomena which at a first glance cannot be explained by applying the 
Fermi liquid theory in its simple form. For example, heavy fermions and the high T c cuprates 
offer such phenomena. In the former they are the metamagnetic behavior and the extremely weak 
antiferromagnetism in the metallic phase. In the latter case they are "spin-gap" behavior and 
"spin-charge separation aspect" , and the controversy about the applicability and limitation of the 
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Fermi liquid theory have been arisen. Itinerant-localized duality model, which has been proposed 
as a quantum phenomenology by Kuramoto and Miyake^ looks for the origin of such phenomena 
with the intermediate energy scale in an explicit treatment of the incoherent contribution which is 
renormalized in quasiparticle in the ordinary Fermi liquid theory. It has been applied not only to 
the heavy fermions systems]!© but also the understanding the anomalous properties of high-T c 
cupratesSBi* 

This model is based on the physical picture that, through the renormalization group evolution 
to the fixed point of such strong correlated electron systems, there is an intermediate stage where 
both the itinerant degree of freedom and "localized-spin" degree of freedom become apparent and 
interact with each other. The physical ground of such a picture is attributed to the fact that the one- 
particle spectral weight has triple peak structure in strongly correlated fermion systems in general 
and large amount of the spectral weight is transferred to the high energy incoherent part while that 
of the coherent one decreases with the absolute value at the Fermi level unchanged. The triple peak 
structure consist of two broad peaks corresponding to the upper and lower Hubbard bands© and 
the narrow quasiparticle peak of Fermi liquid accessible by Gutzwiller's type description. 0) It has 
been explicitly recognized by Kawabata0 two decade ago and demonstrated by recent numerical 
works with d=co technique, d being the space dimension. 00' @@ 0). 

The incoherent part has influence on the physical quantities of the systems even in the low energy 
region where the Fermi liquid description is valid. In particular, the spin fluctuation spectrum is 
affected considerably by the incoherent part, or "localized spins", through a diffusive motion of 
clusters of "localized spins" with short range antifferomagnetic order. This is a salient nature of 
the duality model distinguished from a simple Fermi liquid description where the incoherent part 
is considered to be a simple featureless back ground. A possible importance of the incoherent 
parts was recognized by Ruckenstein and Varma@ in the argument for a microscopic basis of the 
marginal Fermi liquid. We here retain the incoherent part explicitly as a "localized spin" degree 
of freedom which is the irrelevant operator that should be renormalized into the vertex part of the 
itinerant fermion at the fixed point. 

In this paper, we present a microscopic justification of the duality model action for the Hubbard 
model with strong correlation on the basis of the renormalized perturbation expansion of the 
Luttinger-Ward skeleton representation for the thermodynamic potential. In this process, we make 
explicit use of the triple peak structure of the one-particle spectral weight. The organization of this 
paper is as follows. In §2 we develop a microscopic justification of the duality model action and 
in §3 we investigate the role of the high energy process on the fluctuations of the "localized spins" 
which have commensurate correlation through a superexchange-type interaction. Conclusion and 
discussion are given in §4. 
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§2. Justification of duality model - separation of itinerant and localized degree of 
freedom - 

The effective action of the duality model for low energy physics of strongly correlated electron 
systems is written as S 

A = A { + A s + A int , (2.1) 

where 

Af = -EE/t(- ie -)(^i(ien))/ JCT (ien), (2.2a) 

i,j.a n 

As = 2 ^ S i(~' lu rn)(Xoij(^m))Sj{iu m ), (2.2b) 

i,j,m 

A- mt = -^oEE fL(~ ie n ~ iVm)fi(3(ien)Val3 • Sj(il/ m ), (2.2c) 

ia/3 run 

where Af and A s represent the parts of fermions and localized spin consisting of incoherent part of 
fermions respectively, and A{ nt is the interaction term of the fermion with the localized spin. The 
spin degree of freedom, Sj(f m ), mainly comes from the particle-hole exicitations with high energy 
part of fermions, corresponding to the Hubbard bands. However, energies iu m of Sj are restricted 
essentially within the low-energy region comparable to the renormalized bandwidth of itinerant 
fermions. This is assumed by the structure of the dynamical susceptibility Xoijfom) as the energies 
ie n of fermions are restricted by the structure of Gy(ie n ). It is also noted that Sj is a boson field 
describing the quantum fluctuation of the "localized spins" which has been almost quenched by 
high energy process such as the Kondo effect in Ce-based heavy fermions. 

What we want to show below is that the action (2.1) is actually appropriate form for strongly 
correlated electron systems in the low energy region. If the renormalization-flow line of the original 
Hamiltonian of strongly correlated electron systems passes through the point where the duality 
model action can apply in the vicinity of the fixed point, then this duality action is a correct 
effective action of the original one. So we try to show that in the vicinity of the fixed point the 
renormalization drives the original Hamiltonian, which has strong on-site Coulomb interaction, 
toward the form of the duality model action. 

A strategy we take here is to show that the thermodynamic potential in the duality model and 
that of the original Hamiltonian actually coincide with each other term by term under a certain 
condition satisfied by strongly correlated electron systems. The duality model is the effective model 
in the region where the renormalization has been carried out appreciably but not completely. So 
we take the Luttinger-Ward functional form which is the functional of the fully renormalized 
Green function and the self energy of the original Hamiltonian and expand it about the point 
corresponding to the duality model. The thermodynamic potential Y of the original Hamiltonian 
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of the Luttinger-WardllrP form is 

Y = ei6n5 i ln G °(Pi ie «) " *n)G a {p, ie n )] + $[G ff (p, ie„)], (2.3) 

where G a (p, ie n ) and £ CT (p, ie n ) is fully renormalized Green function and self energy respectively 
and &[G a (p, ie n )] is the Luttinger-Ward functinal defined by 

^G a ]=Tj2T,Y,^ n a (p,ie n )G a (p,ie n ), (2.4) 

e„ p,cr n 

where T, n a (p, ie n ) is the ra-th order contribution to the self energy from all the possible skeleton 
diagrams with the fully renormalized Green functions G a (p,ie n ). To avoid an overcounting of the 
same diagram composed of In Green function, the factor *r appears in (2.4). 

We show the expansion of the Luttinger-Ward functional [G a ] up to the 4th order term of the 
Hubbard model in Fig.[l], the solid line is the renormalized Green function and the dotted line 
is the bare Coulomb interaction U. In Fig. [2], on the other hand, we also show the expansion 
of thermodynamics potential 0, of the duality model up to 4th order of coupling constant Ao- In 
Fig. [2], the solid line represents the Green function G a of the fermion of the duality model and 
the wavy line the propageter xo of the "localized spin" , and f2o is the noninteracting part of the 
thermodynamical potential which is written as 

^o = r^^e i5e "lnG (T (p,ie n ). (2.5) 

If the expansion (2.3) of Y, the thermodynamic potential of the original Hamiltonian, is reduced 
to the expansion of f2, that of the duality model, under a certain reasonable condition, we get the 
correctness of the duality model action (2.1) ~ (2.4). To show the equivalence of the expansion Q 
to Y , the diagrams of Y need to be reduced to those of f2 with correct coefficients of corresponding 
diagrams. The expansion of the thermodynamic potential Y contains the fully renormalized Green 
function G a which describes the physics around the fixed point. On the other hand, the expansion 
of f2 contains the Green function of duality model G a which is only partially renormalized. The 
fully renormalized Green function G is expanded as follow. 

G a (p,e) = G a (p,e)+SG (7 (p,e), (2.6) 

where G rT (p,e) is the partially renormalized Green function appearing in the duality model action 
(2.2a) and 5G a (p, e) is the correction due to the renormalization by the interaction (2.2c). Here we 
set up following physical assumptions. 

First we assume that in the energy region, where the duality model can apply, the coherent 
fermions interact only through the spin fluctuations because the charge fluctuations with high en- 
ergy scale are suppressed in strongly correlated electron systems. The diagrams of 6G a (p, ie n ) are 
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given as shown in Fig. [3] where the dotted line represents the partially renormalized Green func- 
tion G a (p,ie n ) and the wavy line the spin fluctuations propagator xo and the dot the fermion-spin 
vertex part Ao- 

Another and the most crucial assumption is that the one-particle spectral weight given by Go- 
already has the triple peak structure in the energy region in question. Namely, the quasi-particle 
peak is narrowed and considerable part of the spectral weight are transferred to the high energy 
region of upper and lower Hubbard band broad peak. The triple peak structure is a characteristic 
feature of the strongly correlated electron systems. Due to this structure we can explicitly dis- 
tinguish the high energy processes of fermion excitations which belong to the incoherent parts of 
the upper and lower Hubbard bands and the low energy processes of fermions which belong to the 
coherent part of quasiparticles. The information that the electron system is in strongly correlated 
region is condensed to the explicit separation of high and low energy part of fermion excitations. 
We can assume this a posteriori. 

Hereafter we consider the diagrammatic structures of the Green function G a (p,e) and 8G a (p, e) 
in detail. The partially renormalized Green function G a (p,e) is written as follow 

G~ 1 a {p, ie n ) = G^oaip, ie n ) - t a (p, ie n ), (2.7) 

where Go a (p, ie n ) is the bare Green function and E CT (p, ie ra ) is a partially renormalized selfenergy 
which is composed only of high energy process. Here the one particle-spectrum, expressed by 
p{x) = —^lmG a (p,x), has triple peak structure already and the low energy coherent parts are 
narrowed while the high energy incoherent parts are broadened. We can separate the high energy 
incoherent process and the low energy coherent process in the spectral representation as 

G a (p,ie n ) = fdx-^)- (2.8a) 
J x-\e n 

dx<^ + / dx*^ (2.8b) 

= G H CT (p,ie n ) + G L CT (p,ie n ). (2.8c) 

The pu{x) is the incoherent part of the one particle spectral weight and the p\ J {x) is the coherent 
one. The Green function G u a (p, \e n ) and the G h a (p, ie n ) are defined as 

G H CT (p,i6„) = fdx^^, (2.9a) 

J X 

G L a (p,ie n ) = [dx-^- (2.9b) 

(2.9c) 



where z is the renormalization factor of itinerant part of duality model, and is given as z = 
(1 — dTj{x) I dx)^ 1 and the £ p is a partially renormalized dispersion. We note that the z becomes 
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already small comparable to the order of z, the renormalization factor at the fixed point. 

Now let us discuss a character of the spin fluctuations which are included in 5G a (p, ie). Due to the 
above assumptions, the spin fluctuations, composed of the incoherent part G H <j(p, ie), dominate in 
the energy region where the duality model can apply. The incoherent parts of fermion are included 
in the irrelevant operator in general. If our standing point is the fixed point of the renormalization, 
the irrelevant operators of the high energy process are renormalized into the vertex parts. However 
once we leave the fixed point, the irrelevant operators will appear no matter how small they are, and 
in strongly correlated electron systems the spin fluctuation processes will be the largest among such 
irrelevant operators. Considering the fact that the particle-hole excitation with high energy fermions 
can carry the low energy, the spin fluctuation process from the incoherent part play an appreciable 
role for the interaction between fermions due to the broad upper and the lower Hubbard bands. 
From the Ward-identity argument^ about the strength of spin-fermion coupling for the impurity 
Anderson model, we can assume that fermions strongly interact with the "localized-spin" . So the 
diagrams of 5G in Fig. [3] can be regarded as we explicitly take out the irrelevant spin- fluctuation 
processes from the selfenergy and neglect other fermion processes owing to their relative smallness 
compared to the spin fluctuations. 

By substituting the Green function G, (2.6), the thermodynamic potential Y, (2.3), can be 
rearranged as follows: 

Y = lnG-£G + $[G] (2.10a) 

= In G (2.10b) 

+ ln (l + ^) (2-lOc) 

- (£ + <5£)(G + 5G) (2.10d) 

+ ^] + Ja\G=Gj'G + j^- 2 \ G= Gj'Gi5'G 2 + ..., (2.10e) 

where 5'G = Gl + 5G, and the summation with respect to frequency, momentum and spin are 
abbreviated. The contribution 5T, represents the excess renormalization on the basis of the duality 
model in which only E has been taken into account. 
The term (2.10b) can be decomposed as 

~,x— A/, n f dx , x i ImG(x) f dx , x ■> ImG(x) 

rVlnGie n =/ — th— tan" 1 ^r 1 + — th— tan" 1 J—^, 2.11) 

^ 1 n) A. GH 2vr 2T ReG(x) JxeL 2vr 2T ReG(x) V ' 

where the symbol x € L and x£H mean that the region of integration is restricted to the low and 
high energy region, respectively. The second term of (2.11) includes only the low energy process 
and is equivalent to lnGL which corresponds to the thermodynamic potential O, (2.5), the non- 
interacting part of the duality model, if G is identified with Gl • 
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The term (2.10-c) is further expanded as 

n=l L 

where we have replaced G by Gl in the right hand side because of the same reason as will be 
discussed below (2.14). The first and the second term of (2.12) are shown in Fig. [4]. 
We next rearrange the term (2.10-d) as 

(£ + <5E)(G + 5G) = £G H + t(5'G) + 5£G H + <5£G L + <5£5G. (2.13) 

The third term of (2.13) is the convolution of <5£ with the low energy components and with the 
Gh high energy ones, and is estimated as 

r£<5£(i en )G H (ie n )= / ^ / ^ ^^^ [/(x)-/(xO], (2.14) 

where /(x) is the Fermi distribution function. The integration value x belongs to the coherent 
low energy region and x' to the incoherent high energy region. Owing to the triple peak structure 
of one particle spectrum, the denominator x — x' is roughly estimated about half of the on-site 
Coulomb energy U in the main contribution of the integral, so this term is negligible. The first 
term of (2.13) gives a constant in the sense of the duality model because it comes from the high 
energy processes. 
Finally, we consider the term (2.10-e). Here the equality 

holds because of the property of the functional $[G], so that the second term of (2.10e) cancels 
with the second term of (2.13). The remaining terms of (2.13) are shown in Fig. [5]. The remaining 
terms of (2.10e) are obtained by replacing Gh's in the skeleton diagrams of <J>[Gf/] with Gl or 
5G by the number of the order of functional derivative of <3?[G] in all possible ways. On replacing 
Gh of the $[Gh] with Gl or 5G, we select the diagrams which can contain the spin fluctuations 
made of Gh, representing high energy process, and discard the diagrams from which we cannot 
extract the spin fluctuation process as irrelevant ones in the energy region where the duality model 
can be applied. For example, of the diagrams shown in Fig.[l] of $[G] , c and f which contain 
only the singlet pair fluctuations are regarded as irrelevant ones and will be discarded. We assign 
the high energy incoherent part Gh for the spin fluctuation part of the diagrams of 3>[G] and the 
spin-fermion vertex parts which correspond to the coupling Ao of the duality model. The remaining 
part of the diagram consists of the low energy coherent part Gl or 5G. 

We illustrate the decomposition of the skeleton diagrams of <$[G] up to the 4th order in U follow- 
ing the above mentioned procedure in Fig. [6], where the thick solid line and thin solid one represent 
Gh, the high energy incoherent part, and Gl, the low energy coherent part, respectively, and the 
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dotted line represents the bare Coulomb potential U, and the wavy line represents the spin fluctua- 
tion propagator xo which represents the longitudinal part Xo or the transverse part Xo~ = Xo+Xq- 
In Fig. [7] we show the spin fluctuation propagator Xo anc ^ Xo~ in terms of Gh up to the second 
order in U. We take up to the first correction of Gl by xo- This is consistent with retaining the 
terms up to the 4th order in U for <3?[G]. The coefficients are attached to the diagrams in Fig. [6] 
because they are important when we compare with those corresponding to the expansion in the 
duality model. It is noted that the bubble of high energy component, in the second diagram of 
Fig. [6-a] , includes both the charge susceptibility and the longitudinal spin susceptibility Xo m this 
level of approximation. However, since the charge fluctuations are irrelevant here, the coefficient 
of Xo should be assigned so as to maintain the isotropy in the spin space by comparing with the 
diagram containing the transverse susceptibility Xq~ > the first diagram of Fig. [6-a] for instance. 

Here, as an example, we explain the diagrams derived from the SG f G \ G=Gu S'G 1 5'G 2 of the 
diagram g in Fig [6]. The diagrams (g-1) and (g-2) are included in the diagram A in the expansion 
of the duality model with the correct coefficient 1/2. Indeed, the diagram (g-1) can be regarded as 
a correction to the first term of Fig. [6-a] which contains the first term of xo given in Fig. [7] . The 
diagram (g-2) contains the two equivalent diagrams which give corrections of two vertices of the 
first term of Fig. [6-a] because the diagram A in the duality model has the spin-fermion coupling Ao 
at both edges of the spin fluctuation xo- The diagrams (g-4) and (g-5) contain the first correction 
to those of (g-1) and (g-2) by the spin fluctuations included in 5G. The diagram (g-5) has two 
equivalent diagrams due to the same reason as for the diagram (g-2). 

Thus the diagrams, derived from the diagrams of Fig.[l-a ~ h] for $>[G] by the procedure 
Sq Sq \ g=Gh S'Gi6'G2, can be rearranged into the diagrams A,B,C and D of Fig. [2] in the dual- 
ity model with the coefficient 1/2, 1, 1/4 and 1/4, respectively. Half of the term corresponding 
to the diagram B cancels with the diagrams of Fig. [4] and Fig. [5]. Indeed, when we sum up the 
diagrams in Fig. [4] ~ [5], corresponding to the terms (2.12) and (2.13) without the second term, 
up to the 4th order in the coupling Ao, we are left with the diagram the same as B of the duality 
model with a coefficient —1/2. So after all we can reproduce the expansion of the thermodynamic 
potential in the duality model, shown in Fig. [2] up to the 4th order perturbation in the spin-fermion 
coupling. This kind of comparison can be extended to higher order perturbation without difficulty 
in principle. It is left for future investigations to show that the above correspondence between 
the duality model and the strongly correlated original Hamiltonian hold to all orders of skeleton 
expansion of 3>[G] in the intermediate stage of the renormalization. 

In conclusion of this section, we have argued on the basis of the renormalized perturbation ap- 
proach that the thermodynamic potential of the Hubbard model with strong correlation can be 
reduced to that of the itinerant-localized duality model. Although we have adopted the Hubbard 
model for conciseness of presentations, the argument can be extended without difficulty to other 
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models such as the periodic Anderson model. 

§3. Spin-fluctuation propagator of "localized spin" exhibiting superexchange-type 
correlation 

In this section we discuss a structure of the spin-fluctuation propagator xo{q,u) on the basis 
of an explicit calculation due to the 2nd order pertubation solution of the two-dimensional (2d) 
Hubbard model. The Hubbard Hamiltonian is written as 

H = J2 e kc{ a c ka + c l- q A'+ql Ck 'i Ck ^ ( 3 - 1 ) 

fccr k,k' q(^0) 

where the conventional notations are used and N is the number of site. The explicit form of the 
one-particle energy is given by means of the hopping term — 1/2 between the nearest-neighbor 
sites as 

6^ = — t(cosk x a + coskyd), (3-2) 

where a is the lattice constant. It is noted that the full bandwidth Wis equal to 41 
The spin susceptibility Xo(<3S i^m) arising from the high energy processes is decomposed formally 

as 

Xo -1 (<?,i^m) = X _1 (i^m) - J(q,iu m ), (3.3) 

where xfam) represents the one-site effect and J(q, ioo m ) intersite effect among the "localized-spin" 
degrees of freedom. If Xo(q,iw m ) has a peak at the wavevector Q = (ir/a,iv/a), it means that the 
nearest neighbor interaction is dominant. By investigating a [/-dependence of the peak height of 
Xo{Q: 0), we can obtain a physical picture for the interaction between the "localized spins". If J(Q) 
scales as 1/U in the limit U /W 3> 1, we can conclude that it can be regarded as the superexchange 
interaction. 

In deriving the duality model in the previous section we have set an important assumption that 
the one particle spectral weight has the triple peak structure. Therefore we have to calculate 
Xo(q, i^m), due to the high energy incoherent part, by means of the Green functions which give the 
triple peak structure to the one particle spectrum. It is well known that the 2nd order perturbation 
with respect to U reproduces the triple peaks structure. So we use the 2nd order perturbed Green 
function G^ 2 \p, ie„), which includes the 2nd order perturbed selfenergy £( 2 )(p, ie„) as an approx- 
imate Green function leading to the duality model. We further approximate xo as the first term 

(2) 

shown in Fig. 7. Namely, we take only the bubble diagram composed of G n (p,ie n ), the incoherent 
part of G^(p, ie n ). The resultant structure of xo will be compared with the approximate form of 
the full susceptibility Xs calculated in the same way but with the use of G^ which contains the 
low energy contributions as well. 
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In order to take out precisely the high energy part of the Green function, we perform the calcu- 
lation on the real frequency axis. The retarded selfenergy as a function of the coordinate R and 
the time t in the 2nd order perturbation with respect to U is given as© 

E^(R,t) = -U 2 [a 2 (R,t)b(R,t) + a(R,t)b 2 {R,t)}, (3.4) 
i 

where the function a(R,t) and b(R,t) are defined as 

«(^*) = ^E ei(P ^ ) /(U, (3.5a) 
p 

b ( R ^ = ^E^'^t 1 - /(&)]. ( 3 - 5b ) 

p 

where £ p = e p — [/,, fi being the chemical potential, and is the Fermi function. The function 

a(R, t) and b(R, t) can be calculated quickly by the fast Fourier transformation (FFT) algorithm. 
The selfenergy in the energy-momentum representation is given by 

POO 

E(p,w)=/ dte^Y e- ip - R rt 2 \R,t), (3.6) 

J ° R 

where both the time integration and the momentum summation are calculated by FFT. The calcu- 
lations are performed with 128x128 sites. We introduce the upper cut-off A for the time integral 
and perform FFT in the interval [0, A] with equally time spaced N e = 1024 points. The cut-off 
A restricts the energy interval to the region [—D,D], where D = irN e /A, and we set A so as to 
D = 5W = 20t 

The chemical potential p, is determined through the relation 



n 



±Y,J dep (2) ( P ,e)f(e), (3.7) 



where n is the electron number density per site and (p, e) is the spectral weight defined as 
/9 (2) (p,e) = --lmG (2 \p, e) 

IT 

1 Im£( 2 )(p,e) 
~ ~7r[e — p — ReS( 2 )(p, e)] 2 + [ImS( 2 )(p, e)] 2 " (3 " 8) 

In Fig [8a] and [8b] we show the calculated results of the one-particle spectral weight of the 
systems with n = 1.0 and n = 0.7, respectively, for various values of U/W. When the Coulomb 
interaction U is comparable with the bandwidth W, the one particle spectrum gets the lower and 
upper Hubbard bands and becomes the triple peak structure, and exhibits the narrowing of the 
coherent peak with increasing U. 

We first calculate Xs(q^), since its calculation is simpler than that of Xoil^) where the low 
energy processes should be excluded. The approximate form of Xs(q,&) given by bubble diagram 
is calculated by means of (p, e) as 

roo 

X .(q,u) =iJ2 dte^ t+ V[X* k+q (t)Y k (t) - Y* k+q (t)X k (t)}. (3.9) 
k Jo 
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where 5 is a positive infinitesimal and the function Xk(t) and Yfc(i) are defined as 

Xk (t) = J dee itf pW(k,e), (3.10) 

Y k (t) = J deeV^M/OO- (3.11) 

The expression (3.9) is the convolution form of Xk(t) and Yfc(i) in the momentum space so that 
it can be calculated quickly by FFT. The calculation is performed at low enough temperature 
T=0Mt=0.0025W. 

The momentum dependence of the approximate full susceptibility Xs(p,u) for a; = is shown in 
Fig[9a], [9b] and [9c] for various fillings. The incommensurate peaks can be seen for all value of 
U/W in the case of filling n=0.7 (Fig[9c]), while they lose their height with increasing the ratio of 
U/W. In the case ra=0.8, the incommensurate peaks disappear even for U/W=1.0. In these figures 
the incommensurate peaks appear in rather large doping for large U/W. This is the effect of the 
approximation that we take but there is a physical reason that the coherent peak which gives the 
incommensurate peak in the susceptibility is narrowed and its weight shifts to the high energy 
incoherent part which is the localized component, so the incommensurate peaks tend to vanish in 
strong Coulomb potential U. We also see that the value of susceptibility decreases with increasing 
U. We explain these below with paying attention to the role of incoherent lower and upper Hub 
bard bands. 

Next we calculate the susceptibility Xo(q, i^m) which is determined only by Gh's representing the 
high energy processes. So, we need an explicit form of Gh- For constructing the Green function 
Gh, we introduce the cut-off which separates the low energy coherent peak and the high energy 
Hubbard bands. Of course the one particle spectrum of the Green function must be the triple peak 
structure for getting such cut-off. 

Here we give Gh as 

G H (p,i£n) = J dx PB ' 

where pn (p, x) is 



X I6 n 

(3.12) 



Ph(p,x) = p(p,x)g(x,C!,c 2 ), (3.13) 

where p(p,x) is usual one-particle spectral function, — ^ImG(p, x), and the function g(x,ci,C2) is 
the cut-off function which cuts the weight of the low energy coherent part that locating in the 
energy region [ci,C2] as 

g(x, d, c 2 ) = 1 + er(a _ e 1 i) + 1 - e r (x „! 2) + 1 , (3-14) 

where the T is a cut-off parameter and we set T = 5t. In order to decide the values c\ and C2, 
we select the saddle points between the the Hubbard peak and the coherent peak. In Fig[10] 
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we show the pu{x) = J2 p Pk(p,x) for the half filled (a) and the n = 0.7 (b) case together with 
P(x) = J2 P p(P,x). 
We give xo{q-,^m) the susceptibility of the incoherent part as 

Xo(<?,iw m ) = -T^2^2G H (k + q,ie n + iuj m )G B (k, ie n ) (3.15) 

£n k 

= ^J dx J dx iu m + x'-x - /(X )] - (3 - 16) 

Non-zero contributions to Xo(dS i^m) arise from the integration regions of x and x' where x belongs 
to the upper Hubbard band and x' to the lower one, or vice versa, because the temperatures in 
question are low enough compared to the energy cut-off separating the coherent and incoherent 
parts. Then we can roughly estimate the denominator iuj m + x' — x in (3.16) for low temperature 
as 1/U and the difference of two Fermi distributions as ±1. This seems to give t 2 /U expected 
for the superexchange interaction. Unfortunately we cannot get the explicit 1/U scaling for the 
peak height of Xo(?>0) in our approximation, although Xo(Q>0) is a decreasing function of U in 
consistent with a general trend expected from the expression (3.16). 

In Fig. [11a] we show the susceptibility of the incoherent parts xo(<2SO), (3.16), for n=0.7 
and U/W=1.0. It is noted that the peak of of Xo(<7jO) is located at commensurate position 
Q = (ir/a,Tr/a). This means that the spin degrees of freedom from the incoherent parts inter- 
act most strongly between nearest neighbor sites. The same point of view about the superexchange 
interaction was suggested by Ohkawa et a/© on the basis of the 1/d expansion and rather rough es- 
timate of the incoherent parts. In our approximation scheme it is apparent that the super exchange- 
type interaction is derived through the electron-hole pair excitations between the lower and upper 
Hubbard bands. The spin susceptibility XoiQ^) appearing in the duality model contains this in- 
formation. 

In Fig.[llb], we show the g-dependence of the full susceptibility Xs(q, 0) calculated by (3.9). It is 
noted that the peak of Xs(q,0) is located at the incommensurate position. This is due to the effect 
of coherent motion of quasiparticle component which reflects the shape of the Fermi surface. It 
is a general characteristic of the duality description that the susceptibility Xo(q), composed of the 
incoherent component, has commensurate correlation while the full susceptibility Xs(q), taking into 
account the coupling with the fermions, has the incommensurate correlation in general. Indeed the 
full dynamical spin fluctuation propagator in the duality model is formally written as 

Xs' l {q,'^m) = Xo _1 (g,iw m ) - A 2 n(g,io; m ), (3.17) 

where the polarization Ii(q, ito m ) of the fermion component is given by 

n(g,iu; m ) = -T^2G L (q + p,iuj m + ie n )G L (p, ie n ). (3.18) 
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Then Xs(Q, ' luJ m) exhibits peaks at incomensurate positions in general, reflecting the g-dependence 
of the polarization U(q,iuj m ). 

It is noted that the incommensurate correlation of Xb(q>w) disappears as the frequency u> exceed 
the cut-off, separating the coherent and incoherent components, leaving only the commensurate 
correlation due to xo(<Z> w )- Such a behavior has been recently discovered in by neutron scattering 
m Cr@ and will be discussed in detail elsewhere. 

§4. Conclusion and discussion 

In this paper we have examined the duality picture diagramatically and the origin of the su- 
perexchange interaction as electron hole pairs excitation between lower and upper Hubbard bands. 
The central assumption about the justification of the duality model is that in the intermediate 
stage of the renormalization, where the duality model can be applied, the one particle spectrum 
already have the triple peak structure and the spin fluctuation process from the incoherent parts is 
relatively larger than other processes, for example than the charge fluctuation. Due to the explicit 
separation between the incoherent/localized and coherent /itinerant parts we obtain naturally the 
"localized-spin" degree of freedoms with neglecting the other fluctuation freedoms. Here a question 
arises whether we can actually separate the coherent and incoherent parts so explicitly. In the 
impurity model like Anderson model and in the d = oo lattice model, d being the space dimension, 
which is equivalent to the impurity model under some conditions the coherent peak of the one 
particle spectrum is narrowed exponentially and the lower and upper Hubbard bands is developed 
explicitly with increasing the on site Coulomb U, we can safely set the above assumption. In the 
d = 2 and d = 3 lattice systems case like Hubbard model, the separation of the coherent and 
incoherent parts is not so clear as in the impurity cases. But we believe the duality model can 
apply to these cases as an effective action when the electron correlation is so strong. 

In §2 we treat the fixed point just like as Fermi liquid but the same procedure is applicable for the 
non-Fermi liquid fixed point if the assumptions we mentioned are satisfied and there no singularity 
appears from the diagrams which we neglected. Actually we do not know the correct analytic form 
of the partially renormalized Green function G(p, ie n ) and the susceptibility of the localized spin 
Xo(q,iu> m ), so we must input appropriate form of G(p, ie n ) and xoil^) which suitably describe 
the physical situation we deal. For example in the reference 5) ~ 7) the weak antiferromagnatism 
of heavy fermion and the normal state properties of cupurate superconductor are successfully ex- 
plained with the suitable form of G(p, ie n ) which include the nesting property. But by considering 
the structures of the incoherent part of fermions which are hidden, for example, in the Fermi Liquid 
theory we get the another properties and explanation of the strongly correlated fermion systems. 
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Figure Captions 

Fig. [1] 

Diagrammatic expansion of the Luttinger-Ward functional 3>[G] up to the 4th order in the 
on-site Coulomb interaction. The solid line represents the fully renormalized Green function G 
and the dotted line the bare Coulomb interaction U. Summation of the momentum, frequency 
and spin index are abbreviated. 

Fig. [2] 

Diagrammatic expansion of the thermodynamic potential of the duality model. The solid line 
represents the partially renormalized Green function G, and the wavy line the spin fluctuation 
propageter xo, and the dot the spin-fermion coupling Ao- 

Fig. [3] 

Diagrammatic expression of SG a (p, ie n ). The wavy line represents xo and the dot-dashed line 
represents Gl, the dot is the spin-fermion vertex part which corresponds Ao of the duality 
model. 



Fig. [4] 

Diagrammatic expression of (6G/G) and (5G/G) 2 
Fig. [5] 

Diagramatic expression of (<5XGl + SUSG) in (2.13). 



Fig. [6] 

Decomposition of the skeleton diagrams for the thermodynamic potential 3>[G] into the di- 
agrams with high and low energy component. The thin line represents Gl, the low energy 
part of the partially renormalized Green function and the thick line Gh, the high energy one, 
and the dotted line the bare Coulomb interaction U. The symbol £ means that the diagrams 
are included in the corresponding diagram (A ~ D) of the duality model shown at the left side. 

Fig. [7] 

Diagrammatic expression of Xoil, 1 ^) and xt'il^)- The solid line represents Gh, the high 
energy part of the partially renormalized Green function, and the dotted line the bare Coulomb 
interaction U. 
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Fig. [8] 

One-particle spectral weight X)pP^(Pi e), given by (3.8), of the 2d Hubbard model for various 
values of U/W with the filling (a) n=1.0 and (b) n=0.7. 



Fig. [9] 

Wavevector dependence of the approximate full susceptibility Xs{q, 0) f° r various values of 
The filling are (a) n=1.0, (b) n=0.8, and (c) n=0.7, respectively. 



Fig. [10] 

One particle spectral weight pn(u), the incoherent part which is indicated by the solid line, 
and p( 2 \co) = ^ p p( 2 \p,w), given by (3.8), which is indicated by the dotted line, for the case 
U/W=1.0. The fillings are (a) n=1.0 and (b) n=0.7 



Fig. [11] 

The wavevector dependence of (a) the susceptibility Xo(q, 0) of the incoherent part and (b) the 
full susceptibility Xs(q,0) for the case U/W =1.0 and n=0.7 



17 



(a) 





U.4U - 




U.OO 




0.30 - 






D" 

w 


0.25 - 








0.20 - 




0.15 - 




0.10 - 




(0,0) 





0.40 






0.35 




o 


0.30 








& 

w 


0.25 




X 








0.20 






0.15 






0.10 





n=0.7 



U/W=0.5 

U/W=0.75 

U/W=1.0 





0.40 - 




0.35 - 




0.30 - 






& 

'w 


0.25 - 








0.20 - 




0.15 - 




0.10 = 



(0,0) 



(0,0) 



(7U,0) 



(0,0) 



(b) 



n=0.8 

U/W=0.5 
U/W=0.75 
U/W=1.0 



j i i i i , I 

M) (jcjc) (0,0) 



